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INFLUENCE  OF  VELOCITY  SHEAR  ON  THE  RAYLEIGH-TAYLOR  INSTABILITY 

I.  Introduction 

An  important  instability  in  a  variety  of  geophysical  phenomena  is 
the  Rayleigh -Taylor  instability.  The  instability  is  an  interchange 
instability  and  is  driven  primarily  by  an  opposing  density  gradient  and 
gravitational  force  (e.g.,  a  heavy  fluid  supported  by  a  light  fluid!.  In 
a  magnetized  plasma,  the  mode  exists  in  both  the  colliBionless  and 
collisional  regimes.  However,  in  many  regions  of  interest,  such  as  the 
ionosphere,  the  plasma  can  also  contain  inhomogeneous  velocity  flows 
transverse  to  the  magnetic  field.  In  fact,  in  the  absence  of  other 
effects  (such  as  gravity,  density  gradients  and  collisions),  this  sheared 
flow  can  give  rise  to  a  transverse  Kelvin-Helmholtz  instability 
(Mlkhallovskii,  1974).  The  purpose  of  this  letter  is  to  investigate  the 
influence  of  velocity  shear  on  the  Rayleigh-Taylor  instability  (Prazln, 
1958;  Chandresekhar ,  1961). 

We  find  that  velocity  shear  can  have  a  dramatic  effect  on  the 
Rayleigh-Taylor  instability.  Namely,  for  a  sufficiently  strong  velocity 
shear,  the  growth  rates  of  the  most  unstable  modes  (i.e.,  those  such  that 
kT,  >  1  where  L  is  the  scale  length  of  the  inhomogeneity)  are 
substantially  reduced;  leading  to  maximum  growth  in  the  regime  kL  <  1. 
Thus,  velocity  shear  has  the  effect  of  preferentially  exciting  a  long 
wavelength  mode  of  the  Rayleigh-Taylor  instability  which  Is  in  sharp 
contrast  to  the  behavior  of  the  mode  in  the  absence  of  velocity  shear. 
This  result  may  explain  the  long  wavelength  oscilllations  (i.e.,  several 
hundred  kms)  of  the  bottomside  F  layer  during  equatorial  spread  F 
(Tsunoda,  1981a;  Tsunoda  and  White,  I9R1;  Kelley  et  al,,  1981)  and  the 
early  time  structuring  of  Injected  barium  clouds  in  the  ionosphere 
(Simons  et  al.,  1980;  Wescott  et  al. ,  1980). 

The  scheme  of  this  letter  Is  as  follows.  In  the  next  section  we 
present  the  basic,  equation  which  describes  the  influence  of  velocity 
shear  on  the  Rayleigh-Taylor  instability.  In  Section  III,  the  various 
instabilities  are  discussed  (i.e.,  Rayleigh-Taylor,  Kelvin-Helmholtz  and 
the  generalized  Rayleigh-Taylor).  Finally,  in  the  last  section,  we 
discuss  the  application  of  our  results  to  ionospheric  phenomena,  i.e., 
equatorial  spread  F  and  plasma  cloud  striations. 

Manuscript  submitted  October  28,  1981. 
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II,  Theory 


The  geometry  of  the  plasma  and  field  configuration  used  in  the 

A 

analysis  are  as  follows.  We  consider  B  •  R  e  , 

A  A  0  8 

g  ■  g  e  ,  n  ■  n  (x)  and  V  *  -  V  (x)  e  .  We  also  include  ion-neutral 
a  x  o  —  o  y 

collisions  (v  -  )  in  the  theory  so  that  both  collisionless  and 

collisional  instabilities  are  considered. 

Hie  basic  assumptions  used  in  the  analysis  are  as  follows.  We 
assume  perturbed  quantities  vary  as  5p  ~  5p(x)  exp(iky-iwt) 

with  a)  <<  and  kr^  «  1  where  0^  is  the  ion  gyrofrequency  and  rTj  is 

the  mean  ion  larmor  radius.  We  neglect  perturbations  along  the  magnetic 
field  (kj  ■  0)  so  that  only  two  dimensional  mode  structure  in  the  xy 
ane  is  obtained.  We  assume  <<  which  is  consistent  with 

ionospheric  F  region  conditions.  Ion  inertial  effects  are  retained  in 

the  analysis,  but  electron  inertial  effects  are  neglected.  We  assume 
quasi-neutrality  so  that  ne  “  *  n. 

A  key  feature  of  our  analysis  is  that  a  nonlocal  theory  is 
developed.  That  is,  the  mode  structure  of  the  potential  in  the  x- 
direction  is  determined  by  a  differential  equation  rather  than  an 
algebraic  equation  obtained  by  Fourier  analysis.  This  technique  allows 
one  to  study  modes  which  have  wavelengths  larger  than  the  scale  size  of 
the  inhomogeneities  (i.e.,  kL  <  1  where  T,  represents  the  scale  length  of 
the  boundary  layer).  In  fact,  this  i3  crucial  to  describing  the  Kelvin- 
Helmholtz  instability  produced  by  a  transverse  velocity  shear 
(Mikhailovskii,  19T4). 

Based  upon  the  assumptions  discussed  above,  the  fundamental 
equations  used  in  the  analysis  are 

3n  /3t  +  7  •  (n  V  )  -  0  (1) 

a  a-a 

V  x 

0  «  F  +  — -  (2) 

—  c 

e  V  x  B 

m^/H  +  1^7)^  -  eK  +  —  g  -  +  m^  O) 

where  a  denotes  species  (e,l:  electron,  ion).  We  substitute  F  ■  -  7^, 
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III.  Results 


To  highlight  the  influence  of  velocity  shear  on  the  Kayleigh-Taylor 
instability  we  first  consider  two  Uniting  cases:  (1)  the  Raylalgh- 
Taylor  Instability  with  no  velocity  shear  and  (2)  the  Kfelvin-Helraholta 
instability  with  no  collisions t  gravity  or  density  gradient.  Following 
this  we  discuss  a  general  solution  of  Eq.  (6). 


A.  Rayleigh-Taylor  instability 

We  take  Vq  -  0  so  that  Eq.  (b)  becomes 


a**  3 Inn  3$  „ 

— *■  + - 2.  _  -  k2 

3xz  3x  3x 


r.~* 

I  0)  + 


,  3  Inn 
1  o 


i\>.  *»  - 

in  3x 


] 


0  (9) 


Equation  (9)  can  be  easily  solved  in  the  local  approximation  (i.e.,  k2L2 


»  k^L2  »  1  where  L  »  (3£nn  /3x) 


-1 

x-x 


'  we  have 


assumed  $(x)  «  exp(ik  x)).  We  obtain  the  dispersion  equation 


w  +  iv.  ta  -  g/L 
in 


0 


OO) 


which  has  the  solution  (Hudson  and  Kennel,  1 975 ;  Haerendel,  1974) 

iv. 


1 v ,  .  i, 

M  - - —  [li(l-  Ag/Lv2  )^ 

-  in 


(11) 


Instability  can  occur  when  g/L  <  0  (i.e..  g  and  L  are  oppositely 
directed).  The  collisionless  and  collisional  solutions  are 


«  »  i  (g/L)1/2 


d)  "  -i  (g/L)  V 


Vin  <<f  <4*/L)  2  02) 


-1 

in 


)1/2  (i: 

^2 


Vin  >:>  (4«/L)  03) 


We  now  solve  Eq.  (9)  numerically  for  a  density  profile 
2  2 

n  ■  n  exp(-x  /’2L  )  +  An.  The  results  are  shown  in  Kig,  1 
which  is  a  plot  of  Y/(g/L)'^vs.  kL  for  An/n  »  0.01  and 


\u 

v^n/(g/L)  0.5,  As  expected,  the  growth  rate  maximizes  in  the  regime 

kL  »  1.  However,  note  that  the  maximum  growth  rate  (Y/ig/L)^2 «  1.1)  is 
slightly  larger  than  that  predicted  from  local 
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theory  (y/(g/L)^*i  1.0  from  Eq.  (11)).  TMe  Is  due  to  the  Gauaaian-like 
density  profile  chosen  for  the  numerical  example* 


B.  Kelvin-Helmholtx  Instability 

We  assume  Vq(x)  *  0  but  take  v^,  -  0,  g  ■  0  and  nQ 

(6)  becomes 


k32V  /&x2 
o 

w  -  kV 

o 


1  ♦ 


0 


-  const*  Equation 


(14) 


which  Is  well-known  (Mikhallovakii ,  1974).  Note  that  In  the 
local  approximation 


o)  -  kV  +  k“l(d2V  /9x2)  (15) 

o  o  x-x 

o 

so  that  the  mode  Is  stable.  In  general,  an  Instability  can  occur  only 
2  2 

when  (3  Vq/3x  )  ■  0  where  *i  <  *Q  <  X2  »  and  x^  and  X2  are  the 

boundaries  (Rayleigh's  the  .ea).  As  an  example,  we  assume  an  equllbrlum 

A 

velocity  profile  V  q  -  Vq  tanh(x/L)  e^.  The  solution  to  Eq.  (14)  is 

shown  in  curve  B  of  Fig  1  where  y/(Vq/’L)  vs.  kL  is  plotted.  The 

instability  is  purely  growing  (i.e.,  u»r  »  0)  and  only  occurs  in  the 

regime  0  <  kL  <  1.  Maximum  growth  occurs  at 

kL  *  0.5  with  y  a  0.18  (V  /L)  (Michalke,  1964). 

o  ■ 


C.  Generalized  Rayleigh-Taylor  instability 

We  now  consider  the  general  case  where  both  the  standard  Rayleigh- 
Taylor  and  the  Kelvin-Welmholtz  instabltties  coexist  (Orazin,  1958). 

Since  the  wavelength  regimes  of  these  two  Instabilities  are  distinct  from 
one  another  (i.e.,  Rayleigh-Taylor  favors  kL  >  1  while  Kelvin-Helmholtz 
favors  kL  <  1),  one  might  think  that  these  "modes"  do  not  affect  one 
another  in  the  linear  regime.  However,  this  is  not  the  case  as  is  shown 
by  curve  C  of  Fig.  1.  Here,  we  solve  Eq.  (6)  for  the  following  profiles: 
n  »  nQ  exp  (-x^/2L2)  +  An 

and  V  -  V  tanh[(x-x  )/L]  e  with  v  /(g/L)^*  0.5  . 

~ o  1/  o  o  y  in 

V  /(gL)'*«  1.0,  An/n  ■«  0.01  and  x  /L  ■  -2.0.  Note  that  the  scale 
o  *  o  o 

lengths  of  the  density  and  velocity  gradients  are  assumed  equal.  The 
position  of  the  sheared  velocity  flow  is  chosen  to  coincide  with  the 
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localisation  region  of  the  Rayleigh-Tayler  Instability.  Aa  In  curve  A, 
we  plot  Y/(g/l.)^va.  kL  and  It  la  found  that  the  growth  rate  maximises 
at  Y  *  0.4  (g/L)^  for  kL  »  0.7  over  the  range  considered  (0  <  kL  < 

10).  The  most  dramatic  feature  of  Including  velocity  shear  in  the 
analysis  is  the  strong  reduction  in  the  growth  rate  of  modes  In  the  short 
wavelength  regime  (l.e.,  kL  >  1).  Thus,  the  main  Influence  of  velocity 
shear  on  the  Raylelgh-Taylor  Instability  Is  to  suppress  the  most  unstable 
waves,  those  with  kL  >  1,  and  to  maximise  growth  In  the  long  wavelength 
regime  (kL  <  1). 

This  effect  of  velocity  shear  on  the  Rayletgh-Taylor  Instability  can 

be  easily  seen  from  local  theory.  For  simplicity,  we  consider  Eq,  (6)  in 

the  limit  k^L^  »  k  ^L^  »  1  and  v.  -  0,  The  dispersion  equation 

x  in 

becomes 

„  ,  S2V  Sinn  SV  Stnn 

0,  -  k  ( — + - )  w  -  g - -  0  (If.) 

SxZ  Sx  3x  Sx 

7  2 

where  8'V  /Sx',  SV  /Sx  and  Stnn  /Sx  are  defined  locally  at  some  point  x  » 

°  I  o  ,  o 

x-.  Writing  S  V  /Sx  ■  V  ,  SV  /Sx  **  V  and  Stnn  /Sx  -  1/L,  the  solutton 
O  O  O  O  0  o 

to  Eq,  (lfi)  is 

w  «*  (v"  +  Vo/L)  +  ~  f(v"/k  +  vVkL)2  +  (17) 

»»  • 

In  the  limit  V  +0  and  V  +  0  Eq.  (11)  is  recovered.  However 
„  ,  o  ° 

for  V  ,  V  *  0  the  velocity  shear  term  is  clearly  stabilizing, 
o  o 

Moreover,  the  stabilizing  influence  is  k  dependent  and  we  expect  that  as 
k  increases  the  influence  of  velocity  shear  becomes  weaker. 

Qualitatively  this  result  is  shown  in  Fig.  1  (curve  C);  the  most  strongly 
suppressed  modes  have  kL  a  1  and  growth  increases  as  kL  increases,  albeit 
small.  However,  owing  to  the  crude  approximation  made  in  obtaining  Eq. 
(17),  good  quantitative  agreement  cannot  be  expected. 
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'/(Vb/L) 
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Figure  1 

Plots  of  the  growth  rate  (y)  vs.  kL  under  various  conditions.  The 

2  2 

profiles  used  are:  n  ■  nQexp(-x  /2L  )  +  An  and 

V  «  V  tanh((x  -  x  )/L].  (A)  Haylelgh-Taylor  instability 
o  o 

with  no  velocity  shear.  Here,  y/(r/L)'2  vs .  kL  is  plotted 

and  we  have  assumed  v^/(g/L)^»  0.5  and 

An/n  »  0.01.  (E)  Kelvin-Helmholtz  instability  with 
o 

no  density  gradient,  collisions  and  gravity.  Here, 

Y/(V  /L)  vs.  kL  is  plotted.  (C)  Generalized  Rayleigh- 
o 

Taylor  Instability  including  velocity  shear.  Here 

Y/(r/L)^vs.  kL  is  plotted.  We  have  assumed 

v  /(g/L)^2-  0.5,  V  /( g/L)^2-  1.0,  An/n  -  0.01 
m  o  o 

and  x  /L  »  -2.0, 
o 
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IV*  Discussion 


We  have  Investigated  the  Influence  of  velocity  shear  on  the 
Rayleigh-Taylor  Instability.  Our  analysis  Includes  gravity*  density  and 
velocity  Inhomogeneities  and  ion-neutral  collisions.  In  general*  the 
Rayleigh-Taylor  Instability  Is  most  unstable  in  the  regime  kL  >  1  where  L 
characterises  the  Inhomogeneity  scale  length.  On  the  other  hand,  the 
Kelvin-Helmholtz  Instability,  driven  by  a  sheared  transverse  velocity 
flow.  Is  unstable  for  kL  <  1.  In  the  presence  of  a  sufficiently  strong 
velocity  shear,  the  short  wavelength  spectrum  (kL  >  1)  of  the  Rayleigh- 
Taylor  instability  Is  strongly  suppressed  and  a  maximum  In  growth  occurs 
for  kL  <  1.  Thus,  velocity  shear  may  cause  a  long  wavelength  mode  to  be 
preferentially  excited;  whereas  in  the  absence  of  velocity  shear  the 
dominant  wave  mode  usually  has  a  shorter  wavelength  determined  by  Initial 
conditions  or  nonlinear  processes.  We  now  discuss  two  possible 
applications  of  this  theory  to  ionospheric  phenomena:  equatorial  spread  F 
and  ionospheric  plasma  cloud  striations. 

It  is  believed  that  the  Rayleigh-Taylor  instability  can  play  a  major 
role  In  the  onset  of  equatorial  spread  F  (see  for  example  Ossakox;,  1979; 
Fejer  and  Kelley,  1980),  After  sunset,  the  density  gradient  on  the 
bottomside  of  the  F  layer  steepens  which  initiates  the  Rayleigh-Taylor 
Instability.  However,  there  are  also  observations  of  (1)  velocity  shears 
existing  In  the  F  region  during  spread  F  (Kudeki  et  al»,  1981;  Tsunoda_et_ 
al. ,  1981)  and  (2)  long  wavelength  (i.e.,  several  hundred  k  is) 
oscillations  on  the  bottomside  of  the  F  layer  (Tsunoda  and  White,  19R1; 
Kelley  et  al. ,  1981).  From  our  theory  we  expect  velocity  shear  to 
preferentially  excite  the  Rayleigh-Taylor  instability  at 
kL  *  0.7.  If  we  take  L  *  25  km  for  the  bottomside  of  the  F  layer  then 
the  most  unstable  wavelength  occurs  at  A  -  250  km  which  is  comparable 
with  observations  (^obs',‘  600  km  Kelley  et  al. , 

1981;  A  .  400  km  Tsunoda  and  White,  1981;  A  -  tens  to  a  few 

hundred  kms,  Tsunoda ,  1981b).  Also,  the  magnitude  of  the  velocity  shear 

necessary  for  this  is  v  /L  =*  2  x  10  Hz  which  is  somewhat  larger  than, 

o 

although  comparable  with,  observational  values  (Kudeki  et  al.,  1981) 

_3 

which  are  V^/L  2  x  10  Hx.  Thus,  the  influence  of  velocity  shear  on 
the  Rayleigh-Taylor  instability  may  explain  the  long  wavelength 


•  / 
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oscillations  of  the  hottomslde  of  the  F  layer.  We  mention  that  gravity 
waves  have  also  been  proposed  as  a  mechanism  to  generate  these 
oscillations  (Kelley  et  al.,  1981). 

Several  aspects  of  the  theory  need  further  comment  concerning  this 
application  to  equatorial  spread  F.  First,  the  velocity  shear  profile 
used  in  the  calculation  is  based  on  observational  data  and  not  on  a  self- 
consistent  equilibrium  model.  An  equilibrium  which  provides  the  observed 
shear  flows  is  beyond  the  scope  of  this  paper.  One  possible  mechanism  to 
generate  the  velocity  shear  is  via  a  neutral  wind  in  the  equatorial  F 
region  when  the  F  layer  is  electrically  coupled  to  a  background  K  layer 
away  from  the  equatorial  region  (Zalesak  et  al. ,  1981).  Tt  is 
interesting  to  note  that  when  such  a  coupling  occurs,  the  plume 
strtuctures  are  tilted  away  from  the  vertical  (Zalesak  et  al.,  1981). 
Thus,  the  tilt  of  the  plumes  can  he  regarded  as  a  measure  of  the  coupling 
between  the  F  and  F  regions,  and  therefore,  as  a  measure  of  sheared 
velocity  flows  in  the  F  region.  Observationally ,  the  largest  amplitude, 
long  wavelength  oscillations  occur  when  the  plumes  are  strongly  tilted 
(Kelley  et  al.,  1981).  This  suggests  that  sheared  velocity  flows  may 
play  a  role  in  their  development.  Second,  the  relative  positions  of  the 
density  and  velocity  profiles  play  a  crucial  part  in  the  "interaction"  of 
the  velocity  shear  and  the  Rayleigh-Taylor  instahility.  The  strongest 
effect  of  shear  occurs  when  the  velocity  shear  is  a  maximum  in  the 
localization  region  of  the  Rayleigh-Taylor  mode.  Finally,  collisions  can 
destroy  the  local  maximization  of  the  growth  rate  in  the  long  wavelength 
regime  if  they  are  sufficiently  strong.  This  indicates  that  velocity 
shear  will  be  more  important  at  high  altitudes  (>  400  km)  in  affecting 
the  Rayleigh-Taylor  instahility. 

Artificial  plasma  clouds  (e.g.,  barium  releases)  in  the  ionosphere 
are  subject  to  a  complex  and  dynamic  evolution.  One  of  the  more  notable 
characrr ristics  is  the  striating  of  the  clouds,  i.e.,  "lingers"  formin 
on  one  side  of  the  cloud  (Rosenberg ,  1971;  Pavis  et  al»,  1974).  In  many 
cases  these  striations  can  be  explained  by  the  x  JB  gradient  drift 
instability  (Linson  and  Workman,  1970;  Zabusky  et  al.,  1973;  Scannapleco 
et  al»,  1976).  However,  recent  shaped  charge  releases  develop  striations 
vary  rapidly  (Simons  et  al.,  1980;  Wescott  et  al.,  1Q80)  and  these 
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Initial  atriattons  cannot  be  explained  by  the  B  drift  Instability 
because  of  its  relatively  slow  growth  rate.  Simons  et  al.  (1980)  have 
proposed  a  kinetic  Instability  driven  by  an  ion  ring  distribution  to 
explain  the  prompt  strlationa  In  the  Ruaro  release.  However,  It  is 
unclear  that  a  kinetic  Instability  can  produce  the  large  density 
perturbiv!  ions  necessary  to  explain  the  structuring  of  the  cloud.  An 
alternative  mechanism  has  been  proposed  by  Pllltp  (1971)  and  Fedder 
(1980)  which  is  based  upon  an  Interchange  Instability;  this  instability 
is  similar  to  the  Rayleigh-Taylor  instability,  but  relies  upon  the 
deceleration  of  the  cloud  ( Scholer,  1970)  rather  than  gravitational 
acceleration.  An  inhomogeneous  electric  field  can  exist  at  the  edges  of 
the  cloud  due  to  polarization  charges  which  produce  a  sheared  transverse 
velocity  flow  (Sydora  et  al.,  1981).  Thus,  our  theory  can  be  applied  to 
the  structuring  of  barium  releases  which  are  injected  across  the  magnetic 
field.  If  the  boundary  layer  is  several  hundred  meters  thick  (L  ~  100  - 
300  m)  then  from  kL  a  0.7  we  obtain  wavelengths  X  ~  1  -  3  km  which  are 
consistent  with  observations.  Moreover  such  a  layer  thickness  yields 
substantial  growth  rates  according  to  Fedder 's  model  (y  ~  10  sec  *)» 

This  final  example  is  largely  suggestive  at  this  time,  yet  is 
sufficiently  encouraging  that  further  investigation  is  warranted. 

In  conclusion,  we  have  shown  that  a  sheared  transverse  velocity  flow 
can  have  a  pronounced  effect  on  the  Rayleigh-Taylor  instability.  For  a 
sufficiently  strong  velocity  shear,  the  short  wavelength  spectrum  of  the 
instability  is  suppressed  and  growth  maximizes  at  kh  <  1.0  where  T,  is  the 
scale  length  of  the  inhomogeneity.  This  result  may  explain  the  long 
wavelength  oscillations  of  the  bottomside  of  the  F  layer  during 
equatorial  spread  F  and  the  prompt  structuring  of  injected  barium 
clouds.  We  emphasize  that  this  is  a  preliminary  report  and  that  a  more 
detailed  analysis  (i.e.,  parameteric  variations)  will  be  presented  in  a 
future  paper. 
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HORTON  AFB,  CA  92409 
(MINUTRMAN) 

01CY  ATTN  MNNL  LTC  KENNEDY 
COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE ,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

01CY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01CY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
01CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
01CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 


SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  XPFS  MAJ  B.  STEPHAN 
01CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
OICY  ATTN  NRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WOPLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 


LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.F.  TASCHYK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 

OICY  ATTN  DOC  CON  FOR  J.  MALIK 

OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

OICY  ATTN  DOC  CON  FOR  J.  ZINN 

OICY  ATTN  DOC  CON  FOR  P.  KEATON 

OICY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUOIJEROUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 
OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 
O' CY  ATTN  DOC  CON  FOR  T.  WRIGHT 


OICY  ATTN  SEA  (SPACE  COMM  SYSTEMS)  OICY  ATTN  DOC  CON  FOR  D.  DAHLCREN 


M.  CLAVIN 


OICY  ATTN  DOC  CON  FOR  3141 


OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 
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SANPIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  ''ASSO 

OICY  ATTN  DC  :  CON  FOR  ».  MURPHEY 
OICY  ATTN  T  C  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 
OICY  ATTN  DOC  CON  DR.  TO  SONG 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI,  RM  5G48,  HQ  BLDG 
WASHINGTON,  D.C.  20505 

OICY  ATTN  OSI/PSW  RM  5P  19 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

(ALL  CORRES:  ATTN  SRC  OFFICER  FOR) 
OICY  ATTN  R.  MOORE 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

OICY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

OICY  ATTN  W.  UTLAUT 

OICY  ATTN  D.  CROMBIE 

OICY  ATTN  L.  ERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  ATTN  R.  GRUBB 

OICY  ATTN  AERONOMY  IAB  G.  REID 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  80009 
OICY  ATTN  I.  GARFUNKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  N.  STOCKWRLL 
OICY  ATTN  D.  OLSEN 


BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIR 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ROMEO  A.  DELIBERIS 

CALIFORNIA  AT  SAN  DIEGO,  UNTV  OF 

P.O.  BOX  6049 

SAN  DIEGO,  CA  92106 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D„B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20/34 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 


ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 


•  I 
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GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

01CY  ATTN  M.H.  BORTNER  SPACE  SCI 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
01CY  ATTN  F.  REIBERT 


INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
01CY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01CY  ATTN  TECHNICAL  LIBRARY 


GENERAL  ELECTRIC  TECH  SERVICES  CO., 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  MILLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.O.  BOX  6770 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 

01CY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN!  SECURITY  OFFICER) 
OILY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
OXCY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANLA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHALL  CROSS 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRFS  ATTN  DAN  MCCLELLAND) 
OICY  ATTN  K.  YEH 


INC.  JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAURAL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
316  STATE  STREET  (P.O  DRAWER  OQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 
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LOCKHEED  MISSILES  *  SPACE  00.,  INC. 

3231  HANOVER  STREET 
PALO  ALTO,  CA  94304 
01CY  ATTN  MARTIN  WALT  DEPT  52-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CY  ATTN  RICHARD  G.  JOHNSON  DEPT  52 
01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  &  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  &  ENGR.  CTR. 

4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALE  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OlCY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OlCY  ATTN  DAVID  M.  TOWLE 
OlCY  ATTN  P.  WALDRON 
OlCY  ATTN  L.  LOUGHLIN 
OlCY  ATTN  D.  CLARK 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OlCY  ATTN  N.  HARRIS 
OlCY  ATTN  J,  MOULE 
OlCY  ATTN  GEORGE  MROZ 
OlCY  ATTN  W.  OLSON 
OlCY  ATTN  R.W.  HALPRIN 
OlCY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OlCY  ATTN  P.  FISCHER 
OlCY  ATTN  W.F.  CREVIER 
OlCY  ATTN  STEVEN  L.  GUTSCHE 
OlCY  ATTN  D.  SAPPENFIRLD 
OlCY  ATTN  R.  BOGUSCH 
OlCY  ATTN  R.  HENDRICK 
OlCY  ATTN  RALPH  KILB 
OlCY  ATTN  DAVE  SOWLR 
OlCY  ATTN  F.  FAJEN 
OlCY  ATTN  M.  SCHEIRE 
OlCY  ATTN  CONRAD  L.  LONGMIRE 
OlCY  ATTN  WARREN  A.  SCHLUETER 


MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  KA  01730 
OlCY  ATTN  JOHN  MORGANSTBRN 
OrCY  ATTN  G.  HARDING 
12  OlCY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

VESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
OlCY  ATTN  W.  HALL 
OlCY  ATTN  W.  FOSTER 

PACIPIC-SIERRA  RESEARCH  CORP 
1456  CLOVF.RFIRLD  BLVD. 

SANTA  MONICA,  CA  90404 
OlCY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OlCY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PROTOMKTRICS ,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

OlCY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OlCY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
OlCY  ATTN  FORREST  GILMORE 
OlCY  ATTN  BRYAN  GABBARD 
OlCY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OlCY  ATTN  ROBERT  F.  LELEVIER 
OlCY  ATTN  WILLIAM  J.  KARZAS 
OlCY  ATTN  H.  ORY 
OlCY  ATTN  C.  MACDONALD 
OlCY  ATTN  R.  TURCO 
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HAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 
01CY  ATTN  BARBARA  ADAMS 


STEWART  RADIANCE  LABORATORY 
UTAH  f*ATE  UNIVERSITY 
1  DR  ANGELO  DRIVE 
BEDFORD,  MA  01730 
OICY  ATTN  J.  ULWICK 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 


RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORR,  NY  10023 
01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 
P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VISIDYNK 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 


SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
OICY  ATTN  L.L.  COBB 
01CY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BABON 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE* 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOODMAN  -  CODE  4180 
Dr.  R.  RODRIGUEZ  -  CODE  4187 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENFSHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/ NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBF.LT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 


NATIONAL  TECHNICAL  INFORMATION  CENTRA 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  H.  MULLANEY 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
PR.  P.  BENNETT 
DR.  F.  KLEVANS 

PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 
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SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
OR.  D.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 


UTAH  STATE  UNIVERSITY 
4TR  AND  HTH  STREETS 
LOGAN,  UTAH  S4322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
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UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-IO,  MS-664 
LOS  ALAMOS,  NM  B7545 
M.  PONGRATZ 
D.  SIMONS 
G.  BARASCH 
L.  DUNCAN 
P.  BERNHARDT 

UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  HILLGARD  AVENUE 
I, OS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A. Y.  WONG 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 
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